Flow Buildup in a Mixed Phase of Quark-Gluon Plasma Plus Hadrons 



m 
o 
o 

(N 
C 

(N 



> : 

ON . 

o\ ■ 
o ■ 
\o : 
o . 

m ■ 

o ; 
^ : 
i 

Oh! 
<D . 
^ ■ 



X 



Judith M. Peters and Kevin L. Haglin 

Department of Physics, Astronomy & Engineering Science, 
St. Cloud State University, 720 Fourth Avenue South, St. Cloud, MN 56301 USA 

(Dated: February 7, 2008) 

Transport processes are considered for their abilities to generate flow in a mixed-phase system of 
hadrons plus quark-gluon plasma exhibiting temperature and density gradients. Diffusive mecha- 
nisms are found to be the leading sources. Even after overcoming viscous effects, conductive and 
diffusive forces generate flow velocities for pions in a model using kinetic theory to within fifteen 
percent of light speed. Since such flow is not observed, this could possibly rule out the mixed phase 
in the evolution of high energy heavy ion collisions. 
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The ultrarelativistic nuclear many-body problem is 
challenging to study in the laboratory because system 
sizes are limited to roughly stable nuclei and the excited 
states are populated for about as long as it takes light 
to traverse the system. And yet, there is a wealth of 
physics information to be gleaned, if only there were 
a way to extract it from the collision zones created in 
heavy ion collisions. Experimental programs at the Su- 
per Proton Synchrotron (SPS) and at the Relativistic 
Heavy Ion Collider (RHIC) are rapidly advancing ow- 
ing to continually developing tools for this extraction [I). 
Various signals of hadronic and subhadronic degrees of 
freedom have been proposed and are being widely used. 
For example, strangeness is expected to be enhanced and 
hidden charm to be suppressed at high energy density; 
spacetime correlations for produced particles provide de- 
vices for spatial and temporal measuring; photons and 
dileptons are the best temperature and spectral probes 
one has; and collective expansion as an observable could 
provide a dynamical measure of the equation of state. We 
study mechanisms which generate flow, and focus on the 
mixed phase. We hope to establish the extent to which 
observed flow signals and existence of a mixed phase with 
gradients are inconsistent with each other. 

Flow is clearly visible in momentum spectra at RHIC 
to a level roughly one- half to two-thirds the speed of light 
depending on the species and the degree of centrality of 
the reaction^. Such explosive evolution in the nuclear 
systems could probe salient features of underlying quark 
and gluon degrees of freedom. The inevitable return to 
hadronic matter from highly excited quark matter might 
proceed through a mixed phase. If so, and if the system 
has even modest temperature and density gradients in 
the mixed phase, we show that conductive and diffusive 
processes give rise to an effective mechanism for gener- 
ating significant flow. So significant that they turn out 
to be the leading sources; much stronger than hydrody- 
namic flow from pressure gradients. Even after overcom- 
ing viscous effects, flow velocities from conduction and 
diffusion in a simple model reach within fifteen percent 
of light speed. These forces would superimpose and add 



to hydrodynamic forces, which are responsible for signif- 
icant flow by themselves 0, Q. One could point to this 
as evidence ruling out the mixed phase with a tempera- 
ture gradient in the evolution of such nuclear many-body 
systems. 

We begin with the system in a state very near equi- 
librium at T c ~ 170 MeV. Local kinetic equilibrium is 
assumed but globally we allow a temperature gradient 
dT/dr ~ —1 MeV/fm. This is the simplest assumption, 
although other profiles could be envisaged. We note that 
a temperature gradient is not inconsistent with hydrody- 
namics since a local maximum temperature different by 
as much as 70 MeV from the average temperature has 
been estimated[5|. By this stage in the evolution of the 
nuclear collisions, probably 5-10 fm/c after maximum 
overlap, the system is roughly spherical in shape with a 
radius of d = 10 fm. For definiteness, let us indicate tem- 
peratures are central temperatures, and we here assume 
zero skin thickness. Two quark flavors plus gluons and 
prehadrons (forming hadrons) are assumed to be in lo- 
cal thermal equilibrium. In practice, this corresponds to 
densities of 5.5, 0.30, and 0.08 per fm 3 for massless par- 
tons, pions and protons, respectively. It has been recently 
pointed out that under these conditions thcrmophoresis 
is responsible for generating significant flowQ. 

Description of irreversible phenomena in systems dis- 
placed somewhat from equilibrium are dominated by 
thermodynamic forces and flows. The flows tend to 
smooth out the inhomogeneities introduced by the forces. 
One such example is heat conduction which is responsible 
for a force per unit area on each hadron given by 
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where n is the number density of partons in the quark- 
gluon plasma, (u) is the mean parton velocity, A their 
mean free path and (E) is their average energy. The 
net thermophoretic force (from conduction) acting on the 
hadrons due to parton dynamics is then given by 



F = irR 2 (F/A) 



(2) 



where R is the hadron radius. Using elementary methods, 
one can estimate from here the result of this radially 
outward force in terms of the asymptotic hadron velocity 
expected in a picture relevant for high-energy heavy-ion 
collisions. In 6] Thoma found the final radial velocity for 
pions to be Vf — 0.81 and for protons Vf = 0.40, although 
see for a correction. 

We go beyond the effects of conduction and begin to 
study other transport processes in the presence of the 
same temperature gradient. If there were a temperature 
gradient, then the default scenario is to support a density 
gradient as well. Particles will of course diffuse from 
regions of high density to regions of lower density. Parton 
diffusion, we argue, is therefore responsible for a force per 
unit area 
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where this time the spatial variation in density drives 
the mechanism. The average momentum appears rather 
than energy trivially, since force is proportional to p • 
Finally, the last effect we include is viscosity. Viscous 
effects arc notoriously difficult to include in hydrodynam- 
ics, but here we use elementary kinetic theory argumen- 
tation to write 
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Viscosity effects, like diffusion, involve average momen- 
tum rather than energy, as was the case for conduction 
which involves energy transport. In the massless limit 
the distinction disappears. Notice also that viscosity in 
this approach is 3/10 of the conductive force and that 
the direction of the viscous force opposes conduction and 
diffusion 8]. Instead of using the simplified approach 
taken here, one could compute collision brackets in the 
first Chapman-Enskog approximation for linearizing the 
full parton kinetic theory 9]. However, the elegant results 
could get lost in the details. We also note that viscosity 
of a quark-gluon plasma has been estimated previously 
using transport theory [Toj| as well as kinetic theory 0. 

Now we explore two different pictures: 1) a system of 
massless partons plus hadrons and, 2) a system of par- 
tons with dynamically-generated masses plus hadrons, 
with m q — gT/\/6 for quarks and antiquarks and m g 
— gT/2 (1 + Nf/6) for gluons, where Nf is the number 
of flavors |l3j. The idea is to compare the preliminary 
and rough estimates of conduction by Thoma to the more 
quantitative formalism discussed here, to compare also 
with diffusion, and finally, to explore finite-mass effects 
arising from the medium. 

It is particularly convenient to first consider massless 
partons since the expressions can then be reported very 
compactly. To set the notation and indicate explicitly the 
partonic states considered, let g q = g q = N s x N c x Nf 
= 2x3x2 = 12 and g g = N s x N c = 2x8 = 16 for the 
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FIG. 1: Magnitudes of forces due to a temperature gradient. 
Heat conduction effects (solid curve), diffusion (dashed curve) 
and viscosity effects (dot-dashed curve) are shown separately. 



degeneracies. Noting that X q = \q = 4/9A s , we find 
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Since C(4) = tt 4 /90 Q, and since \ q = [SCpasT]- 1 
(with the Casmir invariant Cf = 4/3) has been com- 
puted in studies of parton relaxation rates using a hard- 
thermal-loop (HTL) approximation^^ , and finally plug- 
ging in the degeneracies, we arrive at 
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For simplicity, we assume spherical hadrons (prehadrons, 
or extended hadron bubbles) of radius R ~ 1 f m and com- 
pute F = ttR 2 (F/A) for the three transport mechanisms. 
The thermodynamic forces turn out to be proportional 
to the square of the temperature times the temperature 
gradient. However, we allow the strong coupling constant 
to run with temperature and approach an asymptotically 
free value in the T — > 00 limit according to 



<*.(T) 



1 + C* In (T/T c 



(7) 



3 



10 



10 



0.10 




0.20 



10 



-1 



E _ 2 

> 
a) 
o 



10" 



0.10 



I ' 
conduction 
diffusion 
viscosity 

(mossless portons) 




0.20 



FIG. 2: Final flow velocities for pions (solid curve) and pro- 
tons (dashed curve) generated by transport processes. 



where C ~ 0.76, a s (T c ) = 0.48, and T c = 170 MeV. The 
final dependence of the forces on temperature is softer 
than T 2 as shown in Fig. The results suggest that 
diffusion seems to be responsible for the dominant force. 

These radial forces act on prehadrons of mass M and 
consequently generate radial flow. The force is con- 
stant so quite simply the final hadron kinetic energy is 
equal to the net force times the distance d, i.e. K — 
(-Fcond + Fdifi + -Fvisc) d (since F vl y < 0). Using relativis- 
tic kinematic formulae, we findUJ 
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The final flow velocities for pions and protons using d = 
10 fm are plotted in Fig. [21 as a function of temperature. 
At T = T c = 170 MeV the values are 0.86 and 0.49, 
respectively. 

Next we explore the first-order effects of the medium 
to the extent that thermal masses could change these 
results. Relevant quantities appear in the form 
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d 3 p [O]f( P0 ) / d 3 pf( Po ), (9) 



where po = y p + m 2 and where equilibrium distribu- 
tions / = 1/ [exp (/3po) i 1] are taken as appropriate for 
Bosons (Fermions) using the — (+) sign. We calculate 
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FIG. 3: Transport force magnitudes from massive quarks and 
gluons. Conduction is shown as the solid curve, diffusion is 
the dashed curve and viscosity is shown in the dot-dashed 
curve. For each curve, the massless parton limit is plotted as 
the dotted curve immediately above. 
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where x = m (3 — m/T, Ki is the modified Bessel function 
of order i, g is the appropriate degeneracy, n is the den- 
sity of quarks or gluons, ( u ) is the average speed and the 
choice of sign in the summands corresponds to + for glu- 
ons (Bosons) and — for quarks (Fermions). The massless 
limits for the above expressions are naturally recovered 
when x — > since 



lim (zKifz)) = 1 

z—>0 

lim (z 2 K 2 (z)) = 2. 



(11) 



In that case (p) — (E) and (u) — 1. Truncation of the 
sums using only £ = 1 terms corresponds to the Boltz- 
mann limit. Temperature derivatives of the quantities in 
Eq. lfTU|) can be readily computed as needed. 

The first observation regarding the thermal masses is 
that the effects on forces shown in Fig. |3| are quite small 
(at the ten percent level). This insensitivity to medium 
effects on the forces and the expected flow is consistent 
with previous findings^]- The effects are slightly more 
apparent at lower temperatures mainly due to the den- 
sity. 
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FIG. 4: Final flow velocities due to thermal-parton transport 
forces for pion (solid curve) and protons (dashed curve). Re- 
sults for massless partons are also shown for each case (dotted 



For a symmetrical presentation of the results, we in- 
clude Fig. 0] which reports the final flow velocities of 
hadrons (pions and protons) due to transport forces from 
massive partons. Results using massless partons obtained 
previously are shown as well in the figure. 

Before concluding, a remark is in order regarding the 
choice of gradient and system size. Since the final hadron 
kinetic energy is equal to the product of the temperature 
gradient times the system size d, the same results would 
emerge if one instead took d = 1 fm, for then dT/dr 
would be —10 MeV/fm. Consequently, the expected flow 
would be unchanged. 

We summarize our study as follows. Transport pro- 
cesses have been studied in a mixed phase of quarks and 
gluons plus hadrons at or near T c . A modest temper- 
ature gradient is responsible for a net radial force on 
the hadrons having magnitude 11 MeV/fm. The lead- 
ing mechanism or source of flow is parton diffusion which 
pushes hadrons radially outward roughly 9 MeV/fm. Vis- 
cosity introduces a drag force of ~ 1 MeV/fm. In a sim- 
ple picture, these forces are responsible for flow velocities 
0.49c for protons and 0.86c for pions. The first-order ef- 
fect of the medium was included, in which case the forces 
dropped by 10-20 percent. Flow velocities were conse- 
quently scaled back slightly as well, although relativistic 
velocity dependences on forces are clearly nonlinear. 

One might suggest that these result argue for the 



nonexistence of the mixed phase, since flow velocities 
greater than 0.6c are not seen in the data at RHIC. 
Instead, the expansion and evolution of the fireball in 
heavy ion collisions could proceed very quickly through 
the phase boundary and smoothly as in say, a second or- 
der (or higher) phase transition from quark-gluon plasma 
to hadronic matter. 
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